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Abstract 

We discuss relation between three different methods to identify topologically distinct short-range valence-bond 
ground states in spin-i two-leg ladders: the string order parameters, the level-crossing method, and the ground- 
state expectation value of the twist operator. For this purpose we reexamine a phase transition in the frustrated spin 
ladder. We propose a proper bosonization of the string order parameters. Then, these three methods are shown to 
be equivalent reflecting the Z2 X Z2 symmetry breaking of the order and the disorder fields of the two-dimensional 
Ising model. 
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Since the experimental discoveries of several spin 
ladder compounds, there has been considerable inter- 
est in properties of various two-leg spin ladders. The 
ground states of these ladders are well described by 
short-range valence bond (dimer) states with an energy 
gap in spin excitations. In ladder systems, a variety of 
topological configurations of the dimers are possible, 
so that we expect phase transitions between topologi- 
cally different valence bond ground states. 

As a typical model which exhibits the topological 
transition, we consider the antiferromagnetic two-leg 
spin S = i; ladder with frustration described by 



H = ' Sij+i + S 2 ,j ■ S 2 ,j+i) 



(1) 
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+ J±Si,j ■ S2,j + J y. (Sl,j ■ S2J+I + Si,j+i ■ S2,j)], 

where S n ,j denotes an S — | operator at site j of the 
n-th chain. It has been pointed out that there appears 
competition between two phases [2,3]: For J± <C J x < 
J, the ground state is similar to that of the S = 1 
Haldane-gap state described by the Affleck-Kennedy- 
Lieb-Tasaki (AKLT) model, while for J x < J±, the 



system has the same property of the conventional spin 
ladder, and the resonating valence bond (RVB) phase 
including the rung dimer state is stabilized. 

First, we analyze the phase transition by using the 
string order parameter (SOP) which is known to detect 
the hidden antiferromagnetic order in the Haldane gap 
phase of the 5 = 1 Heisenberg chain. In ladder systems, 
the AKLT and the RVB phases are identified by the 
following two SOPs[3], 
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where (• • •} means the ground-state expectation value, 
and a = x,y, z. The composite spin operators for p — 
odd, even are defined by 



'S'odd.j — Sxj + S; 



2,J> 



r.o 



2.J+1- 



(3) 



Then the AKLT and the RVB phases are characterized 
by ©odd + 0, 0« ven = 0, and by C« dd = 0, O e " ven + 0, 
respectively. 

The SOPs can be bosonized by using the relation of 
the S = 1/2 operator exp(mSn j) = 2i5^j. In the con- 
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tinuum limit, however, (D* dd and C| ven have the same 
bosonized form[3]: lim ]x _ y] ^ 00 {e iV ^ 4 '+ (x) e- iV ^'''+ M ) 
where <f>± is the symmetric/antisymmetric boson field. 

In order to distinguish two SOPs, we consider the 
boson representations of the composite spin operators, 

Soddj — ► d x <t>+ H sin V2(/> + ■ cos V2</!>_ (4) 
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cos V2<t>+ ■ sin^20_(5) 



Here S^dd j an d •Seven.j are related by the shift of the 
phase fields, 

4>± -> 4>± ± n/V8. (6) 
This relation indicates that two SOPs are identified as 

0* dd ~ lim {sm[V20+(a ; )] S m[V20 + (2/)]>, (7) 
I x ~ y I — 100 

O z cvcn ~ lim (cos[%/2^ + (x)]cos[ % /2^+(y)]}. (8) 
I*— v\— ►<» 

Note that Eq. (6) reverses the sign of the nonlinear 
term of the symmetric mode of the sine-Gordon model. 
Thus it turns out that 0* dd = Of ven gives the critical 
point. Besides, since the phase field is related to the 
order (en) and disorder (^i) parameter of the two di- 
mensional noncritical Ising models as exp(i\/20+) ~ 
+ i<7i<72[9], two SOPs are identified as 0* dd ~ 
(ai) 2 ^) 2 , z vcn ~ <Mi) 2 <M2} 2 . Thus Eq. (6) for 
4>+ field corresponds to the Kramers- Wannier dual- 
ity transformation Oi «-» /i; for two Ising copies, and 
Codd 7^ (Ccvcn 7^ 0) means the breakdown of the 
Zi x symmetry in the order (disorder) fields. 

On the other hand, it has been pointed out that 
a critical point between two non-degenerate ground 
states such as the AKLT and the RVB phases can 
be identified by the level-crossing of excitation spec- 
tra with different parities under antiperiodic boundary 
conditions [8]. According to the conformal field theory, 
operators appearing in the correlation functions Oj are 
related to the excitation spectra AEj through the scal- 
ing dimensions Xj\ 
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where v is the spin wave velocity. The operators 
sin(V20+), cos(\/20+) appearing in Eqs. (7) and (8) 
are nothing but those corresponding to the energy 
spectra that give the level-crossing point. 

Finally, we discuss another order parameter. Re- 
cently, the author with Todo proposed the ground-state 
expectation value of the twist operator as an order pa- 
rameter to characterize valence bond solid states in 
quantum spin chains [5], 
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The asymptotic form of zl is given by zl = (— l) n [l — 
0(1/1/)], where n is the number of valence bonds at 
boundary. A similar relation was also found in two- 
dimensional dimer systems[6]. Therefore, we obtain 
z L > (z L < 0) for the RVB (AKLT) state. It has 
been proposed that zl defined by Eq. (10) is identified 
by the ground-state expectation value of the nonlinear 
term of the sine-Gordon model[4,5], 

z L oc (cos(V8(t>+)) — > (fJ,i) 2 {^2) 2 - (o-i) 2 (o- 2 ) 2 . (11) 

Thus it turns out that zl has information of two SOPs, 
characterizing the Z? x Z? symmetry breaking in the 
order (disorder) fields by alternating its sign. In case 
that 5*o dd j in zl is replaced by S^ vcn j, the sign of zl 
is reversed [7] due to the relation of Eq. (6). 

From the above argument, three different methods to 
identify the topologically distinct valence bond ground 
states have been shown to be equivalent. We have also 
performed numerical analysis for Eq. (2) using the ex- 
act diagonalization, and observed that phase transition 
point determined by these three methods show reason- 
able agreement. These three different methods are also 
applied to the analysis of phase transitions in the S — 
1 spin ladder with bond alternation[10,ll]. 
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